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Summary. 

Closed  form  solutions  are  obtained  for  the  problem  of  a rigid  asperity 
sliding  with  Coulomb  friction  over  a power  law  viocoelastic  halfspace. 

The  dual  Integral  equations  relating  the  unknown  normal  traction  under  the 
contact  interval  (also  unknown)  to  the  unknown  normal  displacement  outside 
the  contact  interval  are  solved  by  first  reducing  the  system  to  a generalized 
Abel  integral  equation  and  then  appealing  to  the  theory  of  Riemann-Hilbert 
boundary  value  problems.  The  physical  quantities  of  interest  (eg.  the 
coefficient  of  sliding  friction)  are  determined  for  the  three  canonical 
indentors:  a parabolic  pur.ch,  a wedge  punch  and  a flat  punch.  It  is 

observed  that  for  certain  power  law  materials,  singularities  in  the  normal 


traction  field  occur  even  for  the  smooth  parabolic  indentor. 


Introduction. 


Sliding  and  rolling  contact  problems  in  linear  viscoelasticity  theory 
have  received  considerable  attention  by  analysts  during  the  last  two  decades; 
for  exeunple,  see  [3],  [A],  [7]  and  the  recent  review  by  Rutanaprakarn  [6]. 

The  primary  objective  of  these  studies  is  jo  predict  the  dependence  of  the 
contact  pressure  distribution,  friction  coefficient  (or  mechanical  dissipation) 
and  contact  area  on  the  applied  normal  force,  speed  and  viscoelastic 
properties  of  the  slider  (unless  rigid)  and  the  other  body  on  which  contact 
is  made.  It  should  be  noted  that  without  adhesion  and  surface  roughness 
the  sliding  contact  problem  is  equivalent  to  Uiat  of  rolling  contact  for 
rigid  circular  cylinders  and  spheres  pressed  against  deformable  media. 

In  none  of  these  studies  is  an  exact  analytical  solution  found,  except 
the  case  in  which  the  material  is  characterized  as  a so-called  standard 
linear  solid  (SLS)  [3].  This  type  of  representation  is  far  too  simple  to 
fit  actual  viscoelastic  data.  Indeed,  for  rubbery  materials,  the  modulus 
varies  with  time  or  frequency  over  typically  ten  to  twenty  decades  (e.g.  [1), 
[7],  [8]);  at  best,  the  SLS  can  be  fit,  approximately,  to  actual  data  over 
two  decades. 

One  very  important  feature  of  actual  viscoelastic  behavior  is  that  it 
often  obeys  a power  law  in  time  or  frequency  over  much  of  its  range  of 
variation.  On  this  basis,  we  believe  that  it  would  be  very  desirable  to 
obtain  analytical  solutions  to  problems  involving  power  law  materials. 

The  behavior  predicted  with  such  a model  would  certainly  provide  a better 
mathematical  description  of  contact  phenomena  than  with  the  SLS,  and  should 
be  useful  in  guiding  and  checking  numerical  solutions  to  problems  in  which 
real  viscoelastic  data  are  used. 


m 


Inasmuch  as  a power  law  does  not  fit  actual  data  over  the  entire  range  of 
variation,  it  would  be  helpful  to  know  if  significant  error  is  introduced 
in  contact  solutions  because  of  this  lack  of  complete  agreement.  Rutanaprakarn 
[6]  deals  with  this  question,  and  one  can  infer  from  this  study  (especially 
Figs.  5 and  7 in  [6])  that  for  a given  sliding  speed,  any  descrepancy  is 
unimportant  outside  of  a so-called  "dominant  frequency  range".  This  range 
is  limited  to  approximately  four  decades  for  power  law  exponents  which  do 
not  exceed  one-half;  such  is  the  case  for  practically  all  viscoelastic 
solids  (see  [1],  [7],  [8]).  Also  given  in  [6]  is  a procedure  for  selecting 
the  power  law  exponent  and  coefficient  from  actual  data  when  the  latter 
values  are  not  accurately  characterized  by  a power  law  over  any  four-decade 
range. 

In  Section  1 we  consider  the  problem  of  a rigid  indentor  sliding, over 
a viscoelastic  half-space  and  exhibit  the  governing  equations.  In  that 
section  the  mixed  boundary  value  problem  is  reduced  to  a set  of  dual 
integral  equations.  Section  2 contains  an  analysis  of  the  dual  integral 
equations  based  upon  solving  an  equivalent  Riemann-Hilbert  boundary  value 
problem.  The  results  of  Section  2 are  applied  in  Section  3 to  the  contact 
problem  for  three  canonical  indentors;  a flat  punch,  a parabolic  punch  and 
a wedge  punch.  Consistent  with  contact  problems  for  a linearly  elastic 
material,  the  appearance  of  corners  produces  singularities  in  the  normal 
traction  field  under  the  indentor  for  our  power  law  viscoelastic  naterial. 
However,  for  certain  power  law  materials  we  find  jump  discontinuities 
and  singularities  occurring  even  for  smooth  indentors.  We  conclude  Section  3 
with  a brief  discussion  of  complication^  attendent  with  multiple  Indentor 
problems  for  viscoelastic  materials  which  arc  not  encountered  in  elasticity.  < 


Section  4 contains  an  interesting  observation  which  suggests  that  the 
moving  asperity  problem  may  be  cast  In  a variational  framework.  This  is 
accomplished  by  demonstrating  that  for  a parabolic  Indentor  the  contact 
Interval,  which  is  unknown  apriori,  may  be  determined  by  minimizing  a 


suitable  functional. 


Section  1.  Formulation  of  the  Problem. 


The  particular  problem  which  will  occupy  our  attention  in  this  investigation 
is  the  steady  translation  (to  the  left,  say,  with  velocity  U)  of  a rigid 
2-'dimensional  asperity  pressed  against  the  2-dimenslonal  viscoelastic 
halfspace,  y ^ 0. 

Our  ultimate  goal  is  a description  of  the  normal  traction  field  on  the 
contact  interval  and  the  vertical  displacement  of  the  halfspace  outside  the 
contact  interval. 

The  mathematical  formulation  is  as  follows. 

Neglecting  inertia,  the  force  balance  equations  are 


3a  3t 

3x  3y  “ 


3o  3x 

^ + 0 
dy  3x 

and  the  boundary  conditions  are 


■'^j.yCx.o.t)  = k0y(x,o,t) 


—00  < X < 


v(x,o,t)  = f(x  + Ut) 


a(t)  < X < b(t) 


Oy(x,o,t)  = 0 


X < a(t),  X > b(t) 


a(t)  = a + Ut 
o 


b(t)  = b + Ut 
o 


O (x,o,t)dx  = A = total  load. 

a(t)  y 


Here  v(x,y,t)  and  u(x,y,t)  are  the  vertical  and  horizontal  displacements, 

0 is  the  normal  stress  on  x = constan't' lines,  o is  the  normal  stress  on 
X y 

y “ constant  lines,  t is  the  shear  stress,  k is  a given  nonnegative  constant 

xy 
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and  £(x)  is  the  shape  of  the  indentor.  Adopting  the  standard  notation 
3u 


&v  9u  9v  A 9n  3v 

E “ •r->  e * y = -r—  + "Krt  ^ for  the  strains  we  write 

X 3x  y 3y  xy  3y  3x’  3x  3y 


the  viscoelastic  stress  strain  laws  as 


o = 

X 


A(t-T)|^  dT  + 2 


ft 


3e 


G(t~T)^  dT 


o = 


y J. 


A(t-T)-~  dT  + 2 


xy 


3t 

^y^ 


t 3c 

t^(t-T)— ^ dT 


3t 


I C(t-T)-^  dT. 


To  employ  Fourier  transforms  it  is  necessary  to  differentiate  all  the 


equations  (save  for  the  equation  Oy(x,o,t)  = 0,  x < a(t),  x > b(t)) 


with  respect  to  x and  then  adopt  the  Gallilean  variable  s = x ■<*  Ut. 
Define  the  Fourier  transform  of  any  function  g(x,y)  as 


0 00 

g(p,y)  = I g(s,y)  ds. 

>^271  L«, 


Since  A(z)  and  G(z)  vanish  for  z < 0 it  follows  that 


and 


raa  \ 

^3i^) 

= ip  A 

/3o  \ 

1 

= ip  A 

\3s  / 

/3t~  \ 

-i-l 

y/lTf 

= ipG 

ip 

(5 

where  A and  G are  constants.  If  we  take  E to  be  Young’s  modulus  and  v to  J 


Q ^ A 4“  G)  E 

be  Poisson's  ratio  then  it  is  well  known  that  — s—  . Now  if  we 


'V* 


take  G • IpG  and  'X  = ipA  and  appeal  to  standard  correspondence  principles 
(see  [7])  it  follows  that 


U 4-  ?•) 


-V 

E 


.(i-V)  • 


'\>2 


For  a wide  variety  of  materials  v is  a constant  (see  [7]).  Therefore 

-4X 

it  suffices  to  assume  merely  that  E(t)  =A»t  ll(t)  (H(t)  is  the  Heavyside 
function)  to  characterize  our  power  law  viscoelastic  material.  It  should 
also  be  remembered  that  once  s = x + Ut  is  chosen  as  a new  independent 
variable  both  A and  G are  functions  of  t/U  essentially  and  hence 

E , U^<^r(l-g)(ip)“"^ 

*/2ir 


Since 


IpI  -p  « 

G Ik  sgn(p) 

a + 2(j) 

(JT  + 2G)  , 

-1 


Aj^(p) 


and 


Oy(s,y)  = {Oy(p,y);sl 


1 

/2jr 


r 


o (p,y)dp 

0 y 

‘‘P  • 


1 

G _ 2 

A + 2G  1 - V 


we  obtain 


o (s,o) 

y 


ips 


a-l. 


. Jfea-o).  n e ‘^iiKip)  A^(p)dp  ^ 
J 


2^a-;;2)  1-^  /i/2  - ^ 

V 1 - V/ 


!gn(p)) 


By  defining 


(1) 

and 


g(s)  = 2(l-v^)ay(s.o)/(u“<^r(l-a)) 


.7  2 -1 
6 = — tan 
ir 


1/2  - V 


% 

1 - 


the  boundary  conditions  for  v(s,o)  and  o^(s,o)  yield  the  following  dual 


Integral  equations 


(2)  ^ ^{A-(p)>  = f*(s)  a < s < b 

JL  O O 

(3)  ^ ^{Aj^(p)  jpl(ip)“  ^/[l-tan(Bii/2)sgn  (p)];s}  •*  g(s) 


where  g(s)  vanishes  for  s < a and  b < s and  is  unknown  for  a < s < b . 

o o o o 


It  should  be  noted  that  the  parameter  3 satisfies  0 j<  3 < 1 and  that  3 “ 0 
corresponds  to  frictionless  sliding.  The  total  load  condition  requires 


g(s)  to  be  Integrable  on  (a^,b^)  so  that  Eq.  (3)  may  be  inverted  to  obtain 


(4)  “ [1  - i tan  (3iT/2)sgn(p)]  IpI”  (ip)  (2it)' 

Substituting  (4)  into  (2)  yields 


f o 


g(s)e^®^  ds 


(5)  f’(s)  » 5^  ^^IpI  ^(ip)°^  ^ ^ {g(y) ;p);s) 


-i  tan(3Ti/2)  S'*  ^(p  ^(ip)°^  ^ 5**  (g(y) ;p};s}. 


The  right  side  of  (5)  may  be  simplified  as  follows: 


^ ^(IpI  ^(ip)^  “5^  {g(y) ;p};s} 


/OO  f • 

1 


(ip)  e P dp  g(y)e  ^^^dy 


(6) 


b 

f 0 ^ 

“ ^ g(y)t*y  IpI 

/ a — <» 

o 

. rOrair f “ sWf.  . „„  f ° c<t>,Jt 

L^s  (t-s)  •’a  (s-t) 


The  Interchange  required  for  line  (6)  is  readily  justified  since  g(t)  is 
assumed  integrable.  In  a similar  way  it  may  be  shown  that 


^ ~^{p"^(ip)^~°‘^{g(y);pl;s) 

= i sinCaTr)----^^^ 

^ ia  (s-t)^"« 


Substituting  (6)  and  (7)  into  (5)  results  in  the  following  integral  equation 
relating  f'(s)  and  g(s); 


f * (s)tt 

r<l-a) 


cos(a+3/2)iT  rs 


(t-s)^”*^  cos(3tt/2) 


["  _£(t 

^ a (s-t 


.1-a  • 


Equation  (8)  is  a generalized  Abel  type  integral  equation.  The  special 
case  a + 3/2  = 1/2  produces  an  elementary  Abel  equation  which  may  be 
readily  inverted  to  yield 


g(t)  f°-£^da 

r(l-a)  dt  J|,  (5_j)Ci 


However,  obtaining  a solution  for  (8)  when  a + 3/2  i‘  1/2  requires  much  more 
sophisticated  techniques.  One  approach  is  to  map  the  problem  onto  the 
interval  (0,“)  and  use  Mellin  transforms.  Computing  the  resulting 
inverse  transforms  requires  summing  certain  infinite  series  of  residues. 


u; 


The  series  obtained  in  this  way  are  expressible  in  terms  of  hypergeometric 
functions. 

A far  more  elegant  closed  form  solution  of  (8)  from  which  all  physical 
quantities  of  interest  can  be  computed  may  be  obtained  by  solving  a certain 
Riemann-Hilbert  boundary  value  problem  that  is  equivalent  to  (8) . This 
technique  has  been  described  by  Gakov  [2]  and  is  developed  in  the  next 


section. 
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Section  2.  Analysis  of  the  Riemann  Boundary  Value  Problem. 

Define  the  function  of  a complex  variable  ?.  analytic  in  the  complex 
plane  cut  along  the  segme.it  the  real  axis  by 


$(z)  “ [ (z-a^) (b^-z) ] 


-a/2 


g(t) 


'a  (t-z) 
o 


1-a 


dt. 


"•ot/2  ct*"X 

The  multivalued  function  [(z-a  )(b  -z)]  (t-z)  is  defined  and 

o o 

analytic  in  a plane  cut  along  bv  some  branch.  The  density  g(t) 

is  assumed  to  have  the  representation 


g(t)  = 


[(t-a^)(b^-t)l 


a 


where  g*(t)  is  Holder  continuous  on  [a  ,b  ].  Define,  for  convenience, 

o o 


K(z)  = [(z-a  )(b  -z)] 
O 0 


a/2 


Then  $(z) 


R(z)  I 


g(t) 


a (t-z) 
o 


1-a 


dt.  Note  that  4>(z)  is  analytic  in  the  plane 


cut  along  [a  ,b  ] and  that 
o o 


(p(z)  = 


HS  Z ->*  00) 


(9) 


<J)(z)  = 0[(z-a  and 

o 


<l>(z)  = 0[(fa  -z)~^^^]  . 
o 


Let  4>“(x)  denote  the  limiting  values  of  <5'(z)  a.s  z approaches  the  cut  from 
above  and  below  respectively.  It  is  readily  shown  (see  Gakov  [2])  that  on 


12 


<}>^  <x) 


and 


<D  (x)  = - 


1 r. 

“cnri 

g(t)_dt 

’ ° g(t)dt 

R(x)  [ 

a^  (x-t) 

X (t-x)^  _ 

b 

1 

P'  g(t)dt  -aiii 

° g(t)dt 

R(x) 

J . „<,l-oi  ^ 

- ■'a  (x-t) 

0 

X (t-x)  J 

From  these  two  relations  it  is  seen  that 


(10) 


r R(t)dt 

-e  (x)  + 4)  (x) 

'a  (x-t)^"^ 
0 

-2aiii 

e - 1 -• 

R(x) 


and 


(11) 


b 

° _8Ct)  .dt  ^ 

^ (x)  - e 4>  (x) 

, ,1-a 

■'x  (t-x) 

-2cnii  , 
e - 1 J 

R(x) 


Substituting  (10)  and  (11)  into  (8)  we  obtain,  finally,  a 
Rlemann  boundary  value  problem  for  the  determination  of  $(2).  That  is, 
it  is  required  to  find  a function  4>(z)  analytic  in  the  complex  plane  cut 
along  the  segment  subject  to  the  asymptotic  conditions  (9)  and 

whose  limiting  values  <J>  (x)  and  4)  (x)  satisfy  on  the  cut  the  equation 


(12) 


<I>  (x)  * G(x)4>  (x)  + gj^(x) 


. ^ iir(oc+3)  j / \ 2irf’(x) 

where  G(x)  - -*e  and  gj(x)  = r(l-a)R(x)  (1-i  tan  (3ti/2)) 


Once  4>(z)  has  been  determined,  g(t)  may  be  obtained  by  first  computing  the 
limits  $ (x)  and  then  solving  either  of  the  Abel  integral  equations  (10) 
or  (11). 

. » 

The  index,  X»  of  the  Riemann  problem  depends  upon  the  class  of  functions 


J3 


In  which  a solution  is  sought.  In  particuJar,  if  the  solution  is  to  be 
bounded  at  both  endpoints,  X is  -1;  if  it  is  allowed  to  be  unbounded  at 
both  endpoints  X is  +1;  while  if  it  is  to  be  bounded  at  one  endpoinc  and 
unbounded  at  the  other  X is  zero. 

Define  the  auxiliary  function  4>^(z)  by 


^^(z)  = (z-a  )^(z-b 

o o 1 


Im(z)  > 0 


Tz-a  “J^rz-b 

^ n ^ ri 


lm(z)  < 0 . 


The  branch  cut  for 


r—a  ly 


consists  of  a Line  segment  from  z to  a 

o o 

rz-b  lY 


connected  to  the  interval  along  the  x-axis  whereas  for 


~z-b  ~ 
o 

z-z 
L oj 


consists  of  a line  segment  from  z to  b connected  to  the  interval 

o o 

tbo,“).  The  exponent  y is  chosen  to  be  if  is  to  be  bounded 

at  a^  and -f — 2 — 1 otherwise.  The  exponent  y'  is  then  chosen  to  be 
~CY+X)»  where  X is  the  index.  The  function  'l’j^(z)  may  be  determined  by 
solving  the  Riemann  problem  with  boundary  equation 

(13)  $t(x)  = (x-z  )X<l.^(x)  + (x-a  ) ^(x-b  ) ^g,(x). 

i o i o o 1 

The  equation  (13)  is  assumed  to  hold  on  a semicircle  in  the  upper  half-plane 

with  [a  ,b  ] as  diameter  and  z in  the  interior.  If  X = 1 the  solution 
00  o 


to  (13)  that  is  0 as  z is  given  by 


(l-a^)-Y(x-b^)“Y  g^(T)dT  ^ 


T - Z 


I 


where  C is  an  arbitrary  constant  and  the  function  Xj^(z),  the  canonical 


solution  of  the  homogeneous  problem  corresponding  to  (13) , is  given  by 


X,  (z)  “ 1 and  X,  (z)  » (z-z  ) 
X X o 


In  (14)  the  Integral  Is  to  be  interpreted  as  a Cauchy  principal  value. 


If  X “ 0 or  -1  then  C = 0,  and  when  X = -1  in  order  to  obtain  a solution 


it  is  necessary  to  impose  the  orthogonality  condition 


Co 

(T-a^)"\T-b^)  ^ gj^(x)dT  = 0. 


Before  proceeding  to  a discussion  of  the  form  of  the  solution  for  the 


possible  values  of  x we  shall  consider  the  question  of  uniqueness  of  the 


solution  for  the  dual  integral  equations.  Since  the  problem  is  linear  this 


is  tantamount  to  assuming  f'(x)  = 0 for  a^  < x < b^*  and  hence  that  gj^(x)  ■ 0. 


We  shall  allow  only  the  physically  reasonable  solutions  Aj^(t)  which 


correspond  to  a finite  total  load,  that  is  for  which  g(x)  is  integrable. 


The  homogeneous  dual  integral  equations  possess  a non-zero  solution  only  if 


the  homogeneous  Rieraann  problem  corresponding  to  (13)  possesses  non-trivial 


solutions.  From  (14)  non-trivial  solutions  exist  if  and  only  if  X “ 1 


and  are  of  the  form 


(J)^(z)  = C,  4>,(z)  = (z-z  ) 
X X o 


where  C is  an  arbitrary  constant.  The  function  4>(z)  is  then  given  by 


<l)‘^(z)  = C(z-a  )"^“^^^^^^(z-b 
o o 


fiz)  = -Ciz-z  )<“+e-l)/2^-«C«+6) 

o o 


t • 


The  pnase  shift  in  the  expression  for  4>  (z)  results  from  the  two  different 

*-*Y  “Y  * + 

branch  cuts  used  for  [z~z  ] and  [z-z  ] . Tlie  limits  4>~(x)  are  easily 


O’  ■ o 

computed  and  whet;  substituted  into  (11)  yield  the  following  Abel  integral 


equations  for  g(x) 


(16) 


° £m^-C'fx-a  .(3+2a-l)/2 

X (t-x) 


(17) 


f = C'(x-a  ) 


~(3+l)/2^j^  _^j(B+2a-l)/2 
o 


where  C*  and  C"  are  constants. 
From  (17)  we  see  that 


S(t)dt  = 


C'*sin  all 


V 


X (t-a 

O O ^ 


a 


(x-t) 


a 


Observe  that  for  x near  a , 

o 


xxr/  x-(2a+3-l)/2, 

g(t)dt  = 0[(x-aQ)  ] 

a 


from  which  it  follows  that 


g(x)  = 0[(x-a^)’^^“'*’^^^^^]  . 


The  function  g(x)  is  integrable  at  a if  and  only  if  0 < 2a+3  < 1.  Similarly 

o 


we  obtain  from  (16)  that  for  x near  b 


g(x)  - 0[(b^-x)^®’^^/2]  . 


The  above  analysis  shows  that  the  dual  Integral  equations  possess  a 
unique  solution  in  the  class  of  functions  lor  which  g(x)  is  integrable 
if  and  only  if  2a40  ^ 1.  Note,  however,  that  the  solution  is  unique  for 
all  values  of  2a+6  if  g(x)  is  required  to  be  square  Integrable. 

In  the  subsequent  sections  we  shall  have  need  of  the  solutions  of  the 
generalized  Abel  equation  (8)  for  X = -1  and  X = 0.  We  list  them  here 
for  future  reference. 

When  X “ ”1  we  obtain  for  4>“(z) 


41 


! I 

■ f 

^ 3 
--  m 

--  I 


. ^ .(l-a-3)/2.  . .(a+3+l)/2  -iri(a+3+l)/2  b 
. iz-a  ; ) e f o 

2 o_ 


-1 


(i  + tan  ■"(ait/2))r(l-a) 


° (bj-T) 


(t  - z) 


. ^ .(l-a-B)/2.  , .(0ri-&fl)/2  -Tii(a+B+l)/2  b 

^ f °(r-a 

(i  + tan"^(BTt/2))r(l-a)c^^^“'^'^^^  % 


(T  - z) 


The  limits  $ (x)  and  4>  (x)  are  readily  computed  using  the  Plemelj  formulas 
and  when  substituted  into  (10)  and  (11)  yield  the  following  Abel  integral 
equations  for  the  determination  of  g(t) : 


(18)  [ — = r(a)(l+(tan  ^(Btt/2)]^)  [f ' (x)sin(iTB/2)-cos(TrB/2)p(x)] 

U (x-t)-^““ 
o 


(19) 


= r(a)(l  + tan~^(BTr/2)]^)‘^''^[sin[(a+B/2)Tr]f’(x)-cos[(a+B/2)Tr]p(x| 


X (t-x) 


where 
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m 

M 


p(x) 


,(e-i)/2, 


(x-a  )<l-«/2(b  i f “ <'-%) (V"’ 

O o ir  J 

o (t-x) 


-(6+2a+lV2jt(,),,.  I 


• i 


Of  course)  the  above  analysis  is  subject  to  the  orthogonality  condition  (15) 
which  now  takes  the  form 


'1 


(20) 


f “ (T-a  )<6-l)/2(b  -,)-<e<-2“+l)/2f.(T,dT  = 0. 

K ° o 


The  calculations  for  X = 0 are  similar  and,  for  a solution  that  is 


iW 


bounded  at  a^  and  unbounded  at  b^,  result  in  the  following  Abel  equations 


(21)  r (l+[tan‘^(BiT/2)]^)"^^^[f(x)8in(Tie/2)-cos>3/2)p-(x)I 

U (x-tr^  ^ 

o 


(22)  [ s(.0i.t  = r(a)  (l+tan“^(3TT/2)]^)"^''^[8in((a+3/2)TT]f'(x)-cost(a+e/2)Tr]p,(x)l 

^x  (t-x)^"“ 


where 


P^(X)  . (x-a„)<l-««(b^-x)‘^^“-»«  i 


’o  (T-a  )«-«/2(b  -T)-(e+2«-l)/2f. 


(T-x) 


In  the  next  section  we  apply  the  results  obtained  above  to  an  investigation 
of  the  behavior  of  the  solutions  to  the  moving  asperity  problem  for  three 
canonical  indentor  shapes. 


w 

m 


e Ki!  BK-a  ''M  &' 


• 
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Section  3.  Examples. 

In  this  section  we  shall  discuss  the  behavior  of  the  solution  to  the 
moving  asperity  problem  for  three  canonical  indentors,  the  flat  punch, 
the  wedge  punch  ana  the  parabolic  punch,  since  each  one  exhibits  Interesting 
features  of  the  power  law  model.  For  the  sake  of  simplicity  we  shall 
assume  frictionless  sliding,  i,e.  0 0.  Allowing  3 to  be  non-zero  introduces 

no  new  behavior  into  the  solutions  and  represents  a straight  forward 
extension  of  the  techniques  needed  to  treat  the  frictionless  case. 

In  the  physical  problem  the  function  f'(x)  is  assumed  to  have  a 
known  functional  form  for  a given  loading  and  Indentor  shape  on  the 
unknown  contact  interval  [a  ,b  ].  It  Is  required  to  determine  a , b and 
Oy(x»0),  or  equivalently,  g(x).  Moreover,  we  shall  also  compute  the 
coefficient  of  friction  given  by 


(23) 


S 4 


f'(x)Oy(x,0)dx 


where  \ is  the  total  load  given  by 


A = 


0 (x,0)dx. 

y 


It  is  an  interesting  feature  of  viscoelastic  material  models  that  there  is 

a resultant  frictional  like  force  impeding  the  frictionless  sliding  of  a 

punch.  It  is  generated  by  the  bulging  of  the  material  ahead  of  the  asperity. 

We  consider  first  the  flat  punch  or  square  cornered  block  sliding  over 

the  material.  This  corresponds  to  f'(x)  = 0 on  the  contact  interval  [a  ,b  ]. 

• > o o 

If  the  contact  interval  is  to  be  non-degenerate,  i.e.  not  just  a singleton 
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point,  then  g(x)  must  be  a non-trivial  solution  of  the  homogeneous 

generalized  Abel  integral  equation  (8).  By  the  analysis  of  the  last 

section  It  is  apparent  that  an  integrable  function  g(x)  satisfying  the 

homogeneous  equation  over  a non-degenerate  Interval  will  exist  if  and  only 

if  0 < a < 1/2.  Hence  for  1/2  < a < 1 we  are  forced  to  conclude  that  a * b , 

— o o* 

that  is,  the  material  is  in  contact  with  the  flat  punch  only  at  the  leading 
edge  and  the  problem  becomes  that  of  a moving  point  load  I (See  Figure  1.) 

To  analyze  the  flat  punch  problem  for  1/2  < a < 1 it  is  simplest  to 
appeal  directly  to  the  dual  integral  equations  (2)  and  (3)  for  8 » 0. 

Without  loss  of  generality  we  may  assume  the  contact  point  to  be  a^  » 0. 
Therefore,  consider 


Si 


E-1 


P 

i' 


t 


m 


I'i 


(24) 


(25) 


^ ^{Aj^(p);sj  ---  f*(s)  s / 0 


^ ^(Aj^(p)  IpI  (Ip)^  ^;s}  = c6(s) 


where  c is  an  arbitrary  constant  and  6(x)  denotes  the  lirac  measure.  From 
(25)  we  obtain 


Aj^(p)  = cjpl  (ip) 


which  when  substituted  into'  (24)  yields 


f’(s)  = 


r(l-g)c 


-cosan  s s > 0 


ia-1 


s > 0 


The  constant  c is  determined  from  knowing  the  total  load.  It  is  interesting 
to  observe  that  when  a = 1/2,  f'(s)  = 0 for  s > 0 and  that  f'(-O)  =+=>». 

. t 

Hence,  the  displacement  is  constant  (which  may  be  taken  to  be  zero)  past 


■i 


m 


20 


the  contact  point  and  las  a vertical  tangent  as  it  approaches  the  contact 
point  from  the  left.  (See  Figure  1.)  For  1/2  < a < 1 the  displacement 
to  the  left  of  the  contact  point  is  similar  to  that  for  a = 1/2,  whereas  to 
the  right  of  the  contact  point  it  is  no  longer  constant.  Rather,  f(s)  = 0(s^ 
as  a approaches  positive  infinity.  (See  Figure  1.) 

The  situation  for  0 < a < 1/2  is  completely  different.  Now  the  contact 
interval  is  the  entire  block  face  rather  than  merely  one  corner.  From  (16) 
we  see  that  g(x)  may  be  detennined  from 


fi(t)dt  . ^...w  ^-l/2,.  ^a-1/2 

“ 5 ' c cot(ait)(x-a  ) (b  -x) 

i— o o 


'x  (t-x) 


Solving  this  Abel  equation  wc  obtain 

c cos  an  d 


b . ^-l/2,.  xtt-1/2  , 

f o (x-a  ) (b  -x)  dx 


g(t) 


dt 


(.:-t) 


a 


w 


(26) 


cr(i-g) 


/a  r(i/2  - a) 


/u  \^/-u  '“1/2,-..  v-a-1/2 

(b  -a  ) (b  -t)  (t-a  ) 
o o o o 


The  constant  c is  easily  determined  from  the  total  load  condition.  (See 
Figure  2) . 

To  investigate  the  displacement  off  the  contact  interval  it  is  useful 
to  notice  that  equation  (8)  is  valid  over  the  entire  real  axis.  Hence  we 
see  that  for  s > b 

o 


I 


(27) 


f ' (s)  = -cos  an 


■ o s(t)dt 
(s-t) 


whereas  for  s < a 


f’(s) 


' ° -Mt^dt 
•'a,  (t--s) 


In  particular 


-O(s^) 


a+  1/2 


-0((s-b  r" 
o 

-0((a  -s)^^^) 
o 


-0(ls|«) 


8 ■*■  +K. 


s ->■  b + 
o 


s ->■  a - 
o 


s ->  -t”  . 


Note  that  the  displacement  has  a vertical  tangent  at  bctrh  the  leading  and 
trailing  edges  of  the  punch.  (See  Figure  1.) 

We  now  consider  the  case  of  a parabolic  asperity*  i.e.  one  for  which 
f’(s)  ■ “8.  If  a = 1/2  the  problem  is  trivial.  We  obtain  for  g(t)  ■' 


D 

Id  f ° f * (s)ds  , 


The  function  g(t)  will  be  continuous  for  t = b only  if  f’(b  ) 

o o 

if  b a 0.  In  that  case  we  have 
o 


= 0,  i.e. 


g(t)  = ^ ^ 


a < t < 0. 
o 


The  point  a^  is  then  easily  determined  from  the  load  equation.  However, 

we  shall  postpone  this  calculation  until  the  general  case  has  been  considered. 

When  0 < a < 1/2  we  seek  a solution  to  the  Riemann  problem  that  is 

bounded  at  both  a and  b . The  index  of  the  problem  is  then  -1  and  the  function 
o o 

g(x)  may  be  obtained  by  solving  either  of  the  Abel  equations  (18)  and  (19) 

. I 

with  3=0.  Before  describing  the  solution  g(x)  it  is  convenient  to  deterroine 
the  endpoints  of  the  contact  interval  [a^,b^]. 


The  orthogonality  condition  (20)  becomes 


1 
^ Si 


idi  = 0. 

O 


From  (29)  it  is  evident  that  a < 0 < b > i.e.  the  material  wraps  around 

o o 

the  Indentor.  (See  Figure  3.)  The  integral  in  (29)  is  easily  evaluated 


and  yields 


L^m/2,1/2  - = 0 


where  L denotes  the  length  of  the  contact  interval  (b  -a  ) and  B(*,*) 

o o 

Is  the  beta  function.  We  therefore  obtain  the  relation 


L(l/2  - g) 
(1-a) 


The  parameter  L will  be  determined  later  from  the  load  equation. 

To  determine  g(x)  it  is  convenient  to  first  evaluate  the  Cauchy 
Integral  defining  p(x).  It  is  straight  forward  to  show  that 


p(x)  = -X  tan  an  - r(l/2)r(l/2  - a)  (b^-x) 


a+1/2 


■{2F^(l/2,a-l/2;a+l/2;(b^-x)/L}  ~ ^ (iT^rs) 


• ,F,  (1/2,  a+l/2;a+3/2;(b  -x)/L)}. 
i J.  o 


The  derivative  of  p(x)  may  now  be  taken  and  yields 


p’(x)  ~ -tan  an  - r(l/2)r(l/2  - a) (b^-x) 


a-1/2 


•{[(x-a^)/L]^^^[b  a/L  + (1/2  - a)]  - F,  (1/2 .a  - l/2;a  + 1/2;  (b  -x)/L)| 

O O Z A ' 0 
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Since  tan(air)f  * (b^-)  = p(b^-),  the  solution  of  (19)  is  given  by 


(32) 


g(t)  = -ain(aiT)r(a) 


-a 


(x-t)  [sinairf"(x)  - cosairp ' (x)  ] dx . 


Using  (31),  (30)  and  the  fact  that  f”(x)  = -1,  the  integral  defining 
g(t)  in  (32)  may  be  evaluated  and  after  some  manipulation  yields 


(33) 


g(t)  = 


D 

r(l-a) 


(b  -T)~^'^^(x-a  )^'^^‘“dT 
o o 


l-<x 


r(l/2)(l/2-g)L 
r(l-a)r(2-a)r(a+l/2) 


b ( v-1/2,.  vg-1/2, 

; o (t-a  ) (b„-T)  dx 

' o o 


(X“t) 


a 


(See  Figure  4.) 

Nov;  assume  1/2  < a < 1.  For  this  case  it  is  necessary  to  seek  a solution 
of  the  Rieraann  problem  that  is  bounded  at  only  one  endpoint  since  the 
orthogonality  condition  (29)  is  satisfied  only  for  a degenerate  interval 
a^  = b^  and  allowing  this  requires  the  material  to  pass  through  the  asperity. 
It  is  easy  to  see  that  we  xnust  require  the  solution  to  be  bounded  at  a^ 
and  hence  unbounded  at  The  function  g(x)  may  therefore  be  determined 

from  either  of  equations  (21)  and  (22)  with  3=0.  Evaluating  p^(x)  we 
obtain 


> • 
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p;(x)  . -tanai.  t (W)O.  - W2^1/2  - « 


IP 


(1-0)  '•  * '’o 


(x-a  )/L 
o 


L 

* +-(1^ 


Since  for  1/2  < a < 1 the  orthogonality  equation  (29)  does  not  hold,  some 

other  condition  must  be  Imposed  to  specify  a^  and  b^.  It  is  natural  to 

choose  the  points  a and  b so  as  to  enhance  the  endpoint  behavior  of  g(x). 
o o 

Clioosing 


(34) 


^ - (^)  ^ 


forces  g(x)  to  be  continuous  at  x - a^  while,  as  with  the  case  0 < a < 1/2, 


choosing 


/1/2-ax 


forces  continuity  at  b^.  However,  the  condition  (34)  is  physically  impossible 

since  it  predicts,  among  other  things,  a total  load  of  the  wrong  sign. 

The  function  g(x)  must  be  positive.  From  equation  (8),  it  is  clear  that 

if  1/2  < a < 1 then  f*(s)  must  be  positive  on  [a^,b^].  Since  f'(s)  = “S  ‘.this 

means  that  b <0,  whereas  line  (34)  results  in  b > 0.  Hence,  we  conclude  that 
o - o 

the  condition  (30)  obtained  for  0 < a < 1/2  is  valid  also  for  1/2  ^ a < 1. 

Note,  then  that  for  1/2  < a < 1 the  material  leaves  the  indentor  before  the 
apex.  (See  Figure  3.) 

Given  relation  (30)  it  is  straight  forward  to  show  that  the  formula  (33) 
obtained  for  g(t)  when  0 < a < 1/2  is  valid  for  1/2  < a < 1.  It  should  be 


observed  that 


« f 
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Hence  it  is  apparent  that  when  a > 1/2  g(s)  is  monotone  decreasing  on 

(a  «b  ).  However,  when  0 < a < 1/2,  g(s)  has  a unique  maximum  occurring 
o o 

at  the  point 


= -(S) 


(See  Figure  4.) 

The  behavior  of  the  displacement  profiles  may  be  deduced  from  equations 
(27)  and  (28) . It  should  be  noted  that  the  displacement  has  a continuously 
varying  tangent  over  the  whole  line. 

Mention  should  be  given  to  the  question  of  uniqueness  for  the  parabolic 
indentor.  As  was  demonstrated  in  the  last  section,  when  1/2  £ a < 1 uniqueness 
is  assured,  whereas,  when  0 < a < 1/2  the  solution  to  (8)  is  non-unique 
in  the  class  of  integrable  functions.  However,  only  one  of  the  solutions 
is  continuous.  Adding  to  (33)  any  non-zero  solution  of  the  homogeneous 
problem  results  in  a displacement  profile  whose  tangent  becomes  vertical 
upon  approaching  the  contact  interval.  This  is  obviously  physically 
impossible  since  it  corresponds  to  the  material  passing  into  the  indentor. 

As  was  indicated  earlier,  a quantity  important  in  predicting  the  response 
of  a viscoelastic  material  tc  loading  is  the  friction  coefficient  given 
by  (23).  For  the  parabolic  asperity  is  p; sily  computed.  In  particular 
if  f'(x)  = -X  on  (a^,b^)  then 


t'j 


i 


9 


^f  * (l-a)(3-a)  ^ • 


For  completeness  we  shall  indicate  wlmt  the  solution  for  the  parabolic 
asperity  looks  like  when  the  friction  «terra  is  included.  For  g(t)  we  obtain 


0 

I 


a 
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f °(b  -X 

1-a)  J ^ 0 


o 


^ n > (h  -v'»^P+2a-l)/2,  ^ x-(B+l)/2 

r((B4-l)/2)(l-2a~B)L^  ° 

2r(l-a)r(2-a)r({l+(^2a)/2)  J a 


The  trailing  contact  point  becomes 


b . W2 

o (1-a) 


where  L comes  from  the  load  condition 


X = U°‘d;r((3-B-2a)/2)r((B+3)/2)L^~^ 

2(l-v^)(l-a)r(3-a) 


Finally  for  we  obtain 


r - (a+g)(2-g)  , 
^f  ■ (3-a)(l-a) 


It  should  be  noted  that  the  parameter  that  governs  the  essential  behavior 
of  the  solution  in  this  case  is  a + B/2  rather  than  a. 

We  next  consider  a wedge  shaped  punch  for  which 


f(x)  = 


X < 0 


X > 0 


where  and  are  positive  constants.  Assume  first  that  a ^ 1/2.  Since 

g(t)  must  be  non-negative,  we  see  from  (8)  that  a and  b must  be  such  that 

o o 

f'(x)  > 0 on  (a  ,b  ).  Hence  b < 0. 

~ o o o ~ 

For  a = 1/2  the  problem  is  again  trivial.  The  function  g(x)  is  given  by 


(See  Figure  6.)  Recall  that 


1 


! 


1 


f’(x) 


Sj  ,”0 


n 


a 


(b  -t)~^''^(t-x)“^^^dt 
o 


t > b 


a < t < b 
o o 


t < a 


Obviously,  from  physical  considerations  we  must  take  b = 0.  The  leading 

o 

contact  point  is  now  easily  determined  from  the  load  equation.  (See  Figure  5.) 
It  should  be  noted  that  the  stress  is  singular  at  the  trailing  contact  point 
(the  apex  of  the  wedge)  and  has  a jump  discontinuity  at  the  leading  point. 

(See  Figure  6.) 

When  1/2  < a < 1 we  seek  a solution  to  the  Riemann  problem  that  is  bounded 

at  a and  unbounded  at  b . As  before  we  solve  the  Abel  equation  (22)  for 
o o 

g(t)  where  3 “ 0 and  Pj^(s)  is  given  by 


p^(s)  = 6 tan  an  - 


5r(l/2)r(l/2-a)(b 

0 

TTl’d-a) 


-F,(l/2,a  - l/2;a  + 1/2; (b  -s)/L). 
i 1 o 


This  reduces  to  solving 


(35) 


b 

f 0 


g(t)dt 


s (t-s) 


1-a 


= h(s) 


where 


(36)  h(s)  - 


S,r(a)r(l/2)(b  ■ “ 

I o 

r(l-a)r(l/2  + a) 


„F,(l/2,a  - l/2;a  + 1/2; (b  -s)/L), 
^ i o 


The  derivative  of  h(s)  is  easily  taken  and  yields 


h'(s)  =- 


S r(a)r(l/2)L^"“(b  -s)“  " ^^^(s-a 
1 o o 

r(l-a)r(a  - 1/2) 


Solving  for  g(t)  we  obtain 


b 

■n  air  [ ° h*  (s)d£ 

J 

•'t  (t-s) 


S,r(a)sln(aTr)r(l/2)  A (b  -s)^  " ^^^(s-a  )'^^^ 
1 , i-ai  o 0 


iTr(l-a)r(a  - 1/2) 


(s-t)' 


(See  Figure  6.)  It  is  easily  seen  that  again  we  must  choose  b^  = 0.  Indeed, 
recall  that 


mi 


If 

.. / X r(l-a)  ° g(t)dt 

f ' (s)  = -cos  aiT  ^ "i 

" ia. 


s > b . 
o 


For  s > b f ’ (s)  is  a decreasing  function.  Since  f'(s)  = S.  for  a < s < b 
o loo 

we  conclude  that  f*(s)  < S,  when  s > b . This  is  possible  only  if  b =0. 

1 o o 

Therefore,  when  1/2  < a < 1 the  material  leaves  the  wedge  indentor  at  the 
apex  in  contrast  to  the  parabolic  indentor  for  \diich  it  leaves  in  front 
of  the  peak.  (See  Figure  5.) 

The  load  equation  is  readily  computed  and  yields 

^ ^ U^cVr(3/2  - g)  ^1-^a 

2irr(2-a)(l-v^)  ■ • 


from  which  L,  and  hence  a , may  be  determined.  Note  also  that  for  1/2  < a < 1 

o — 

the  friction  coefficient,  C^,  is  just  the  slope  of  the  leading  face  of 


the  wedge,  i.e. 


It  should  be  observed  also  that  the  stress  is  singular  at  both  a^  and 
b^.  In  particular 


(0((s-a 

{ -1/2 
(0((b  -s) 


s -*•  a + 
o 

s -►  b - 
o 


(See  Figure  6.) 

Assume  now  that  0 < a < 1/2.  We  may  observe  first  that  in  this  case 
> 0 > a^,  i.e.  the  material  must  xnrap  around  the  indentor.  (See  Figure  5.) 
Indeed,  from  (37)  we  conclude  that  for  0 < a < 1/2,  f*(s)  is  nonnegative 
when  s > b , from  which  it  follows  that  b >0.  Moreover,  if  b =0  then 
we  may  appeal  to  the  analysis  for  the  case  1/2  < a < 1 to  conclude  that 
g(t)  must  solve  (35)  with  h(s)  given  by  (36).  Inverting  (35)  yields 


g(t)  = - 


sin  air  d f ° h(s)dt 


n dt 


it  (s-t)« 


Observe  that  h(s)  =*  0[(b  -s)^  for  s ->•  b - which  implies  that 

o o 

-1/2 

g(t)  * 0[(b  -t)  ] for  s -V  b - . This,  together  with  (37)  shows  that 

o o 

f'(s)  = 0[(s-b^)°^  for  s b^+  , from  which  it  follows  that  the 

displacement  possesses  a vertical  tangent  upon  approaching  b^.  The  last 

statement  is  a physical  impossibility.  Hence,  b >0. 

o 

To  solve  for  g(t)  we  seek  a solution  of  the  Riemann  problem  that  is 

bounded  at  a and  b . It  is  required,  then,  to  solve  equation  (19) 

° . » 

subject  to  the  orthogonality  constraint  (20).  However,  for  the  wedge 


punch  with  a^  < X < b^,  the  solution  of  equation  (19)  is  not  expressible 


i 
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tn  terms  of  elementary  functions  but  may  be  estimated  numerically.  Moreover, 

simple  expressions  for  and  b^  are  not  obtainable,  though  it  is  possible 

to  calculate  a^  and  b^  numerically  as  follows.  Define  a new  parameter  0 to 

be  o * -a  /L.  The  orthogonality  relation  (20)  becomes 
o 


0 » [ 2“^^^(l-z)““  ‘ ^^^dz  - S [ z"^^^(l-z)"“  " ^^^dz. 

■'o  ^0 


The  load  equation  Is 


. uV,r(i-a)  f°  ... 

^ — g(s)ds 

2(l-v  ) a 

o 


From  (18)  with  8 = 0 we  obtain 


'a  (x-t)' 
o 


= - r(a)p(x) 


where  p(t)  is  given  by 


li 


p(0  = H 


Hence, 


^o  (T-a  )‘^^^(b  -T)*^  " 

o o 

a_  (t  - t) 


f ’ (T)dT 


f = f ° 

■*  a l'(l-<»)  ^ a (b  -t)^ 

o o 
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TTr(i>a) 

o 


-T)-“  - f 

O 


'>0  (t-a 

o o 

a (t  - T) 


I 

1 

a-“)  J. 


° ^-l/2,.  ^x-a  - 1/2^, 

T(T-a  ) (b  -T)  dr 


r(l-a)  " o 

a 

o 


r(l-a) 


E 


1/2.,  x-a  - 1/2.  - 1/2,,  .-a  - 1/2  , 

z (1-z)  dz  - z (l-z)  dz 


Line  40  follows  from  an  obvious  change  of  variables  and  repeated  applications 
of  (38).  To  determine  a^  and  b^,  one  need  only  solve  (38)  numerically  for  O 
and  then  use  (39)  and  (40)  to  obtain  L. 

As  was  Indicated  above,  It  Is  not  possible  to  express  g(t)  in  terms 
of  elementary  functions.  Once  a^  and  b^  are  known  g(t)  ma;  oe  approximated 
numerically.  However,  due  to  space  limitations  we  shall  not  address  this 
question  here. 

We  remark  at  this  time  that  our  analysis  yields  as  a special  case, 
namely  a « 0,  the  solution  to  the  problem  of  an  indentor  pressed  into  an 
elastic  half  space.  Indeed,  setting  a = 3 = 0 in  the  above  expressions 
for  g(t)  produces  the  results  obtained  by  Muskhelishvilli  [5]  for  the 
elastic  contact  problem.  It. should  be  noted  that,  as  expected,  the  friction 
coefficient,  C^,  vanishes  when  a = 0 and  the  pressure  distribution  is 
symmetric  on  the  contact  interval  for  a symmetric  punch. 

As  another  remark  we  include  a brief  discussion  of  the  problem  of 
multiple  asperities.  In  particular,  assume  that  several  indentors  are 
pressed  into  the  half-space  and  are  moving  with  the  same  constant  velocity. 
Assume  that  the  indentors  have  shapes  given  by  f ^^(x) , . . . ,f^(x)  and  that  the 


I 


’*  » 
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total  pressure  under  each  of  the  punches  is  known.  Hence,  we  assume  as 
known  the  n parameters 


r 


0.  (x,o)dx 


i “ 1, . . . ,n 
,th 


where  is  the  unknown  contact  interval  for  the  i*""  indentor. 

Letting  g^(t)  denote  the  restriction  of  g(t)  (given  by  (1))  to  the  interval 
it  is  straight  forward  to  show  that  f^(x)  and  g^(x)  i = l,...,n  are 
related  by  the  system  of  dual  integral  equations 


52^“^A^(t);x)  = f^(x) 


X G i * 1, . . . ,n 


where 


ItjCit^'^'^jx]  “ Xj  (x)g 

i=l  ^i  ^ 


(x) 


(1  X e I 
X,  (x)  = { 

i (O  X ^ I 


For  simplicity,  assume  n = 2.  Analogous  to  (8)  we  obtain  the  coupled 
system  of  integral  equations  for  gj^(t),  i = 1,2 


trf^*  M 
Til-a) 


(41) 


r(l-a) 


r 1 8;^(t)dt  rx  gj^(t)dt  r 

— - COS  aiT  + 

••x  (t-x)  (x-t)^  J 


2 g (t)dt 

7r;i-a  * ^1 

a,  (t“x) 


= cos  an 


^ 1 gj^(t)dt 

.■’2 

g2(t)dt 

gj^(t)dt 

<x-t)^"«  ^ i 

x 

- cos  an 

1 

The  elastic  contact  problem  for  multiple  indentors  Is  no  more  difficult 
to  solve  than  for  a single  punch.  The  situation  for  the  viscoelastic  half- 
space is  considerably  different.  For  the  elastic  half-space  it  Is  possible 
to  reduce  the  problem  to  that  of  solving  a single  Rlemann  boundary  value 
problem.  Because  of  the  occurence  of  branch  points  on  each  of  the  Intervals 
no  such  reduction  is  possible  for  the  viscoelastic  problem.  Rather 
one  obtains  a coupled  system  of  Rlemann  boundary  value  problems  fot  which 
no  simple  uncoupling  seems  to  exist.  In  general,  the  system  (Al)  must 
be  solved  numerically. 

There  is,  however,  one  special  case  for  which  (Al)  is  easily  inverted. 

In  particular,  when  a“  1/2  the  second  equation  uncouples  and  is  easily 
Inverted  to  determine  g^it) . With  g2(t)  known,  the  first  equation  may  be 
solved  for  gj^(t).  The  four  parameters  a^^,  b^,  i = 1,  2,  are  still  to  be 
deduced.  The  contact  points  b^  for  the  trailing  indentor  are  determined 
Independently  of  the  first  indentor  since  the  first  indentor  produces  no 
effect  on  the  second.  The  obvious  coupling  in  (Ala)  is  reflected  in  the 
choice  of  a^^  and  b^^  and  the  resulting  relative  displacement  of  the  two 
punches.  For  example,  if  f^(x)  = -x  and  f^(x)  = -(x-c)  for  0 < c then 


and  a^  and  b^^  may  be  calculated  from  the  equations 

b = - I r (c-t)^^^(t-b  )"^^“dt 

1 71  1 

J n 


Section  4.  Observations  on  determining  the  contact  interval. 

In  the  last  section  the  criteria  for  del t rinining  the  contact  interval 
for  each  of  the  examples  depended  upon  the  particular  indentor  shape 
considered.  No  general  principle  other  than  preserving  the  physical 
integrity  of  the  asperity  was  invoked.  For  more  general  indentors  it  may 
be  very  difficult  to  apply  this  criterion  due  to  the  complicated  nature 
of  the  integrals  involved. 

In  attempting  to  circumvent  this  difficulty,  we  conjectured  that  for 
a fixed  loading,  X , the  correct  choice  for  a^  and  b^  should  correspond 
to  a minimum,  or  at  least  stationary,  point  of  when  considered  as  a 
function  of  a^  and  b^.  In  particular,  we  tested  the  conjecture  by  considering 
the  following  constrained  minimization  problem  in  the  case  of  a parabolic 
Indentor  without  friction: 


(42)  minimize  YCj(L,a^) 


ii 


O 


g(s)f'(s)ds  subject  to 


g(s)ds 


Y 


where  g(s)  is  given  as  the  formal  solution  of  (18)  for  0 < a < 1/2  and  of 

'V2 

(21)  for  1/2  < a < 1 and  Y = 

u“6r(l-a) 


The  problem  (42)  may  he  solved  by  a standard  Lagrange  multiplier  argument. 


We  present  the  analysis  for  the  case  1/2  < a < 1;  the  other  cases  may  be 
treated  similarly. 

It  may  be  shown  that  for  f'(s)  = -s 


”i 


T(L,a  ) 
o 


I) 

' o 

g(s) (-s)ds 
a 

o • ' 


;| 
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Conclusion . 

Morland  [4]  has  shovm  that  for  the  relatively  general  case  of  a Kelvin 
material  the  contact  pressure  distribution  over  a smooth  Indentor  vanishes 
at  the  leading  and  trailing  edges  of  the  contact  area,  just  as  it  does  for 
elastic  media.  In  contrast,  however,  we  have  found  that  for  the  parabolic 
asperity  the  pressure  at  the  leading  edge  has  a finite  discontinuity  when 
a ■ 1/2  and  an  infinite  discontinuity  when  1/2  < a < 1.  This  behavior  is 
apparently  closely  related  to  the  fact  that  the  Initial  value  of  the  power 
law  relaxation  modules  1s  Infinite  whereas  the  Initial  value  of  the  Kelvin 
modulus  is  finite. 

Nevertheless,  the  utility  of  our  results,  aside  from  the  fact  that 
closed  form  solutions  are  provided,  Js  renect(>d  by  the  formula  for  the 
friction  coefficient.  This  coefficient,  like  the  stress  intensity  factor 
for  crack  problems,  has  proved  to  be  a useful  characterization  of  the 
response  of  viscoelastic  material  to  sliding  contacts.  When  one  recalls 
that  linear  elasticity  theory  predicts  a physically  unrealistic  singular 
stress  field  at  a crack  tip  but  a finite,  acceptably  accurate  intensity 
factor  the  analogy  becomes  even  more  appropriate. 

The  other  interesting  deviation  from  classical  behavior  that  should 
not  go  unmentioned  is  the  algebraic  displacement  to  ±“  of  the  material 
behind  (ahead  of)  the  aspe'ity  when  a > 1/2.  This  is  to  be  contrasted 
to  the  two-sided  logarithmic  displacement  to  for  an  elastic  half  space. 
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Figure  1.  Displacement  fields  for  a flat  punch  for  a < 1/2 (- 
a = l/2(-*  — ♦);  a > l/2( ). 


Figure  5 


Displacement  fields  for  a wedje  indentor  for  ct 
a » l/2(-*-*-*);  01  > l/2(— — ). 


< 1/2  ( ); 


Figure  6 


Normal  traction  fields  over  the  contact  interval  for  a 

wedge  indentor  for  a < i/2( );  a = l/2(-»~*-.); 

a > l/2( ).  The  interval  is  only  representative, 

it  is  not  intended  to  imply  that  the  same  interval  Is  correct 
for  all  a. 


